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CHARACTERIZATION OF THE CLASS OF 
CANONICAL FORMS FOR SYSTEMS OF LINEAR 

EQUATIONS 


JC NDOGMO 

Abstract. The equivalence group is determined for systems of 
linear ordinary differential equations in both the standard form 
and the normal form. It is then shown that the normal form of 
linear systems reducible by an invertible point transformation to 
the canonical equation = 0 consists of copies of the same 
iterative equation. Other properties of iterative linear systems are 
also derived, as well as the superposition formula for their general 
solution. 


1. Introduction 

One of the major difficulties associated with the study of systems 
of linear ordinary differential equations, and in particular their sym¬ 
metry properties is the large number of parameters that occur in their 
expression as arbitrary coefficients. Therefore, most of the general re¬ 
sults concerning the symmetry properties and point transformations of 
this type of equations are limited to second- and third-order systems 
[Dial [31111 El El [7]. Moreover, the study of these systems has also been 
often limited in the recent literature to the case of constant coefficients 

[SlEl E]- 

In a slightly wider context, Gonzalez-Gascon and Gonzalez-Lopez 
estimated in HD] upper bounds dm,n for the dimension of the symme¬ 
try algebra of systems of ordinary differential equations (odes) of the 
general form 

y(") = F(a;,y(„_i)), y G M'", (1.1) 

where y{n-i) denotes y = y(a;) and all its derivatives up to the order 
n—1. It was thus found in the exact bound value djn ,2 = m^+4m-|-3, 
and the more general value dm,n = + nm+ for n > 2. The value 

of dm ,2 had been obtained earlier by Markus mu in an independent 
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study. All these results for dm,n are in fact extensions of the well-known 
results obtained by Lie himself [12] for di i.e. for scalar equations. 

In the case of second-order or third-order equations of the form fll.ip . 
their group classihcation has been studied for linear equations with ar¬ 
bitrary coefficients mm, and most often with constant coefficients 
[HI inilG] . In the case of arbitrary coefficients, this study has generally 
been reduced in the recent literature to systems of two equations while 
the dimension of the system is often arbitrary for equations with con¬ 
stant coefficients EE]. Another important aspect of equations of the 
form fll.ip is their linearization criteria, and several independent results 
have also been recently obtained about these criteria for second-order 
equations [a Ea n E]. 

A study more directly related to the topic of this paper concerns 
systems of the form fll.ip with the right-hand side equal to zero, and 
which will henceforth be referred to as the canonical form. In the case 
of scalar equations. Lie [15] showed that for n > 2 an nth order ode is 
of maximal symmetry if it is reducible by an invertible point transfor¬ 
mation to the canonical form . More recently, Krause and Michel [TB] 
proved the converse of this result in the case of linear equations and 
showed that in fact a linear scalar equation is of maximal symmetry 
if and only if it is iterative. Some recent results have also been ob¬ 
tained about linear scalar equations of maximal symmetry, concerning 
their transformation properties [min], their coefficient characteriza¬ 
tion [misniin], their first integrals or Hamiltonian forms [2I1122]. 

In this paper we completely derive the equivalence group for linear 
systems of odes in both the standard form and the normal form, where 
normal form here means the standard form in which the coefficient of 
the term of order n — 1 has vanished. We then use this result to show 
that a linear system of odes is reducible by an invertible point trans¬ 
formation to the canonical form if and only if its normal form consists 
of copies of iterative equations with the same source parameters. Sev¬ 
eral other results concerning systems of iterative equations and their 
general solutions are also obtained. The base field for all differential 
equations considered will be assumed to be an arbitrary field K of 
characteristic zero. 

2. Iterations oe linear equations 

We shall denote by A[xi ,... a differential function of the vari¬ 
ables Xi, ... ,Xp. Let T = R^ -|- S be an ordinary differential opera¬ 
tor, where R = (rj^) and S = {sij) are m x m matrices with entries 
the smooth functions rjj and of the independent variable x G K, 
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and where R is invertible. Linear iterative equations are the iterations 
R, s] = = 0 of the hrst-order linear ordinary differential 

equation \['[y] = Ry' + Sy = 0 where y = y(x) G K”, and given by 

vl>’^[y] = [\[/y], for n > 1, where \I>° = / (2.1) 


is the identity operator. This definition is a natural extension to 
vector equations of that given for scalar equations for instance in [121 
in]. A linear iterative equation of a general order n thus has the form 

= Kl y<”> + Kl y<”-‘> + Kl y<”-"> + ■ ■ • + A'r ‘ y' + KZ y = 0, 

( 2 . 2 ) 

where the are mxm matrices. By writing down the recurrence re¬ 
lation among the similar to that found for iterative scalar equations 
in [13 P. 3], the equality = R"^ is immediately found. In the sequel 
1x1 matrices will be identified with scalars while for scalar equations 
we shall usually use the notations r, s, and y instead of R, S, and y cor¬ 
responding to the vector equations counterpart. In the case of scalar 
equations, since the coefficients are scalars and hence commutative, the 
identities 


Ki = r 


n—1 


ns -|- 


Ki=r 


n—2 




3sr' -|- rr" -|- 


3n-5 


.'2 


are readily established. Moreover, the well-known change of the 
dent variable 


(2.3a) 

(2.3b) 

depen- 


y = gw^ where g satishes Kl^ g + n g' = 0 (2.4) 

maps fl2.2p into its normal form. This transformation however simply 
amounts to relating r and s in such a way that = 0. Therefore, for 
given source parameters r and s of the operator T, an nth order scalar 
linear equation in normal form 

yin) ^ ^2 yin-2) + ... + /[i yin-J) + ... A^y = 0 (2.5a) 


is iterative if and only if 


Ai = 


Ki 


(2<J<n), 


Ri=0 

It follows from (12. 3 a^ that the requirement that 
alent to having 


(2.5b) 


0 holds is equiv- 


s = --(n - l)r', 


( 2 . 6 ) 
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and this shows in particnlar that any iterative scalar eqnation in normal 
form can be expressed in terms of the parameter r alone, i.e. depends 
on a single arbitrary fnnction. Moreover, for any fnnction ^ = ^(x) 
which may also be a matrix fnnction, by setting 

^K)(x) = (K'(X)]^ - 2«x)r(x)) (2.7) 

it follows from fl2.3bp that in fl2.5p we have 


A 2 = A{r), and more generally 



( 2 . 8 ) 


In fact, as already noted in [HKITj, the coefficients depend only on 
the fnnction A{r) = A 2 and its derivatives. For simplicity, it will often 
be convenient to denote the coefficient A 2 of the term of third highest 
order in fl2.5ap simply by q. 

If we set r = V? for a certain nonzero scalar fnnction n, the expression 
for A{r) = A{v?) in fl2.7p is mnch simpler and rednces to —u''/u. 
Setting q = A{r) is thns eqnivalent to letting n be a nonzero solntion 
of the eqnation 

l/" + qi/ = 0, (2.9) 

which is referred to as the second-order sonrce eqnation for fl2.5p . 
Differential operators generating linear iterative eqnations in the nor¬ 
mal form fl2.5p with the coefficients = 24^[q] expressed directly as 
a fnnction of q and its derivatives have been obtained in ini and an 
improved version appears in |22] . The most general form of linear iter¬ 
ative scalar eqnations in normal form is given for example for the low 
orders three and fonr by 


0 =2q'?/ + A(\y' + y'" = 0 (2.10a) 

0 =102/'q' + lOqy" + ?>y (Sq^ + q") + y^^\ (2.10b) 

while the corresponding second-order eqnation with the same sonrce 
parameter is given by (12.91) . 


3. Equivalence group 


The gronp of eqnivalence transformations is well known for the gen¬ 
eral linear nth order scalar ode [231123] . For homogeneons eqnations, 
it consists of invertible point transformations of the form 


x = f{z), y = g{z)w, 

(3.1a) 

for the standard form, and 


x = f{z), y = 

(3.1b) 
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for the normal form, where / and g are arbitrary functions and c an 
arbitrary nonzero scalar. For the general linear nth order vector ODE 
which without loss of generality may be assumed to be in the homoge¬ 
neous form 

0 = y(-) + + ■■■ + 3’^ -1- • • • + S" y, (3.2) 

where y G K” and the = {B^-{x)) for n hxed are m x m matrices, 
the corresponding extension of fid.lbll can be guessed. But since such 
a result does not seem to be available in the literature, we completely 
derive it in the next theorem, as well as its normal form counterpart. 

Theorem 1. The equivalence group of the general linear system fl3.2l) 
consists of invertible point transformations of the form 

x = f{z), y = Tw. (3.3a) 

These transformations reduce for the normal form of 03.21) in which 
the coefficient of the term of order n — 1 has vanished to 

x = f{z), y = f{zY^-^'>/‘^Cw. (3.3b) 

Here, T = {Tij{z)) is an mx m matrix and f as well as the entries Tij 
ofT are smooth functions, while C = {Cij) G K” is a constant matrix. 

Proof. Let x = /(2;,w), and y = T(^, w), where w = {wi,... ,Wm), 
/ is a scalar function and T = (Ti,... ,Tm) a vector function, be an 
invertible point transformation for 03.21) . The resulting transformation 
of 03.21) yields a polynomial in the Wj and their derivatives, and the 
coefficients of the corresponding nonlinear monomials must vanish. If 
we denote by p > 1 the highest exponent of wl in this expansion, then 
the coefficient of w'h’ has a term of the form TiBf.^{fyj.y, where for 
a function F = F{ai ,..., a„) we use the notation Fa^ = Fa^aj = 
■TT—a —, etc. Due to the arbitrariness of the coefficient 5L this term 
must vanish, and the invertibility of the transformation then implies 
that fwj = 0, and this holds for all j. Consequently, f = f{z) depends 
on ^ alone. Transforming again the equation with this new expression 
for /, the vanishing of the coefficients of monomials of degree two in 
the polynomial expansion of the transformed equation shows that 

for z, A;,/= 0,..., m. 

Consequently, we must have 

Ti{z, w) = Tij{z)wj + ai{z), for z = 1,..., m, 

for some functions Tij and a* of z. Moreover, from the fact that the 
equation is homogeneous it follows that a* = 0 for all z. This shows that 
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an equivalence transformation of the equation must have the stated 
form fl3.3ap . Conversely, a direct calculation shows that each invertible 
point transformation of the form fl3.3al) leaves the family of equations 
fl3.2p invariant, and thus fl3.3ap represents indeed the a general element 
of the equivalence group of fl3.2p . 

When the general linear system fl3.2l) with = 0 is transformed 
by a group element of the form fl3.3al) . the resulting coefficient of 
is an m X m matrix with entries 




nTlj 


;) T^J^‘ 


for i, j = 1,..., m. 


f/n y/n+l 

Solving each of the equations iff = 0 for Tij shows that an equivalence 
transformation of the normal form counterpart of fl3.2p must be of the 
form fl3.3bp . and a direct calculation shows that each transformation of 
the stated form (I3.3bp leaves indeed the normal form of (13.21) invariant, 
and this completes the proof of the Theorem. 

□ 


It should be noted that expressions for the scalar form fl3.ll) and the 
vector form fl3.3p of equivalence transformations are the same, and in 
particular fl3.ip is just as expected the one-dimensional version of fl3.3p . 


4. The class of canonical forms 

In the study of symmetry properties of systems of linear second-order 
ODES of maximal symmetry, Gonzalez-Lopez showed in a seminal paper 
PP that in its normal form y" + ^(x)y = 0, the coefficient A{x) of y 
must be a scalar matrix. It was also proved in that paper that the 
second-order system of equations is of maximal symmetry if and only 
if its symmetry algebra is isomorphic to s[(m-|- 2, M), assuming M to be 
the base field of the equation. 

It was already known by then mm that the maximal dimension 
of the symmetry algebra for second-order systems of the form fll.ip is 
+ 4m -|- 3 and attained for equations in canonical form y" = 0. For 
linear systems fll.ip of order greater than two, it seems to be known only 
for third-order equations |2] that the maximal dimension is m^-|-3m-|-3 
and is reached only for equations equivalent by point transformations 
to the canonical form y'" = 0. 

The result of [1] about the normal form of second-order systems of 
maximal symmetry was obtained by a direct an intricate analysis of 
the determining equations of the symmetry algebra. By making use 
of equivalence groups, we shall give a similar characterization for the 
class of equations fl3.2p reducible by an invertible point transformation 
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to the canonical form. The class of canonical forms will henceforth be 
often referred to simply as the canonical class. 

Theorem 2. A linear nth order system of the form fl3.2p is reducible 
by an invertible point transformation to the canonical form ^ if 

and only if its normal form has expression 

0 = H-(4.1a) 

where is the identity matrix of order m and the Af = Af{z) are 
scalar functions such that 

0 = + Al + ... + AIW, (4. lb) 

is a scalar linear iterative equation in the dependent variable w = w{z). 

Proof. Given that an equation of the form (13.21) is reducible by a point 
transformation to the canonical form if and only if its normal form can 
be obtained from the canonical form by a transformation of the form 
fl3.3bD . to prove the theorem it suffices to show that the transformed 
version of the canonical form = 0 under fl3.3bH is of the required 
form fl4.1al) . Under fl3.3bD and multiplication from the left by the ma¬ 
trix f'^ 2 ~C~^, the canonical equation = 0 is transformed indeed 
into 


0 = Al Im W^"' H-h A” /„ w, (4.2a) 

where 

0 = + Al{z)w^'^-^^ + • • • + A^{z)w (4.2b) 

is exactly the corresponding transformation of the scalar version y^^'^ = 
0 under fl3.1bp . But since by a result of [16], a scalar linear equation is 
reducible by an invertible point transformation to the canonical form 
if and only if it is iterative, fl4.2bp is iterative and this completes the 
proof of the theorem. □ 

Note that the factors occurring in the right hand side of fl4.1a|) 
may clearly be omitted, and systems of equations of this form are 
termed isotropic. Roughly speaking. Theorem [2] states that the normal 
form of any linear equation in the canonical class consists of copies 
of the same iterative scalar equation. This theorem shows both the 
strength of normal forms for systems of linear equations and the sim¬ 
plicity of the expressions of equations in the canonical class. Indeed, 
even in their normal form not all systems of linear equations fl3.2l) have 





















the simple form (I4.2al) of a completely uncoupled system in which the 
coefficient matrices are not only diagonal, but scalar. Moreover, in its 
most general expression, the normal form of the linear system fl3.2p de¬ 
pends on (n—1) X arbitrary functions, while in fid.lap it depends only 
of n — 1 arbitrary functions, and in fact just on q and its derivatives. 
On the other hand, a system of linear equations reducible to canonical 
form by a point transformation need not be isotropic. Indeed, under 
the equivalence transformation fl3.3ap . the transformed version of the 
generalized free fall equation y” = 0 takes the form 


0 = w" -|- 

f" ' 

2T-^r - 

f 1 

w + 

f" 

rjn — \rp!f J rjn — Irjnf 


[ /' \ 




and thus, in standard form, it is a coupled system that effectively 
depends on eight arbitrary functions, and not just formally on one such 
functions. 

In terms of the iteration operator d' in (12. ip . Theorem [2] says that 
an equation of the general form fl3.2l) belongs to the canonical class 
if and only if its normal form is an iterative equation in which the 
coefficient matrices R = r/^, and S = sIm are scalar matrices and the 
scalar functions r and s are related as in fl2.6l) by s = — (n — l)r'/2. 

Corollary 1. A linear system of the form fl3.2p belongs to the canon¬ 
ical class if and only if its normal form consists of uncoupled iterative 
equations of the form 

^”[2/*] = 0, fori = l,...,m (4.4) 

where 4/ = r{x)-^ -\- s{x) is a differential operator and the scalar func¬ 
tions r and s satisfy s = ^{n — l)r'. 

Proof. This is just a restatement of the above remark according to 
which any equation in the canonical class must be iterative with the 
specihed scalar matrices, and follows by expressing such a system com¬ 
ponentwise. □ 

Although an equation from the canonical class does not consist of un¬ 
coupled systems of identical iterative equations, it can be transformed 
into such a system. 

Corollary 2. A linear system of the form fl3.2p belongs to the canonical 
class if and only if it can be transformed to a system of uncoupled 
iterative equations of the form 

^"[l/*] = 0, fori = 

where 4/ = r{x)-^ -t- s{x) is a differential operator and the source pa¬ 
rameters r and s are arbitrary scalar functions of x. 
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Proof. Indeed, any equation of the general form (13.21) can be trans¬ 
formed into the normal form by a change of the dependent variable of 
the form y = Q(a:)w where the invertible matrix Q satishes 

BlQ + nQ' = 0. (4.5) 


By Corollary [H the equation is in the canonical class if and only if 
this transformed version consists of uncoupled iterative equations of 
the form (14.dh . Then, the additional invertible change of dependent 

variables w = P(x)h, where P = exp 4 B{t)dt 


Im is a scalar ma¬ 


trix and B{t) a scalar function will transform the uncouple system of 
iterative equations with same source parameters to another uncouple 
system of (identical) iterative equations with same source parameters 
in standard form, and this proves the corollary. □ 


It follows from Corollary [2] that iterative equations generated by 
operators of the form 4/ = -|- s with source parameters the scalar 

matrices R and S belong to the canonical class. It is however clear 
that contrary to the scalar case, not all operators T generate equations 
belonging to the canonical class. Indeed, for a general operator 4/ of 
the form 12. Il and y = y(a;) G K” we have 

4>^[y] = Ry" -h (rr' RS -h SR)y' -h (RS' S^)y, (4.6a) 

4/3 [y] = R3y"' 

-h (rr'r -h 2 r^r' -h R^S SR^ RSR)y" 

-h (rr'^ -h R^r" RR's -h 2 r^s' -F rs'r rsr' (4.6b) 
-F RS^ -h S^R srr' -F SRs)y' 

-h (rr's' -F R^s" RS's -F RSS' SRS' s3)y. 


Using the fact that reducing the normalized iterative equation (with 
leading coefficient /„) to normal form amounts to relating R and S in 
such a way that the coefficient of the term of second highest order 
vanishes, it turns out that condition fl2.6p also holds for commuting 
matrices, just because fl2.3l) holds for such matrices. That is, for com¬ 
muting source parameters R and S, reducing the normalized equation 
to normal form amounts to setting 


s = 




(4.7) 


In particular R and r' should commute in such a case. We denote by 


d (n - 1) , 

^'nor = ' R 

dx 2 


(4.8) 


2 
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the corresponding iteration operator. For m = 2, the most general form 
of matrices R commuting with their derivative r' is given by 

where the Xj are arbitrary constants while the a* = ai{x) for z = 1,2 
are some arbitrary functions. When condition (14. 7 j) is satished, the 
coefficient ^(r) of y for the normal form of ('hnor[y] = 0) is clearly 
given by fl2.7l) . However, for a matrix R = (“ ^) of the form fl4.9l) . the 
entries of ^(r) = (aij) are given by 

an = (52 (a'2 + _ 2aa") - + d( - <^7' («' + <i') 

+7 (q;'2 + + 25a") — a (o'y' + -j'S' — 25j ")) 

an = (52 {a'p' + p'5' - 2ap") - 2p^-f5" + /?( - (5 (d'y' + 5'^) 

+7 (a'p' + p'5' + 25p") - a (p'i + 5'^ - 255")) 

(4.10) 

021 = 7 (-^ (a'^ + P'l') + /? (aV + - 27«'0) 

—a'y (a;'2 + /?'y' — 25a" — 2/3y") + a^ {a'^' + ^'5' — 2(5y") 

022 = 7 W + S’) + P {p'i + (5'2 - 2yd")) 

-ay {a'p' + p'5' - 25p" - 2p5") + a^ {p'i + (5'2 - 255 "), 

where a^- = 4ajj(d7 “ o:6)^. It is thus clear that even when the ma¬ 
trices R and S commute and satisfy fl4.7l) . the corresponding matrix 
71 (r) is not a scalar matrix as required for the normal form of the 
corresponding equation to lie in the canonical class. Moreover, the 
nonlinear system of equations ai2 = 0 = 021 representing the necessary 
condition for M(r) to be scalar is not obvious to solve directly. By 
considering some particular cases of this system and by denoting u and 
V arbitrary nonzero solutions of the second-order source equation fl2.9p . 
it turns out that the resulting matrices R satisfying fl4.10l) . i.e. of the 
form fl4.9p and for which M(r) is scalar are given by 


when d = 7 = 0, or 


R = 




(4.11b) 


when the entries of R are all scalar multiples of each other, and where 
the kj are arbitrary scalars. A direct calculation shows that equations 
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generated by iteration operators corresponding to these two matrices 
lie indeed in the canonical class. 

On the other hand, for any choice of q in fl2.9l) and any operator d'nor 
in fl4.8p for which R is a matrix of the form fid.lip , when one applies to 
the generated equation the substitution 

fb) ^ ^ 1-2 ^ fQj. j > 2, and f = u, or v, (4.12) 

where Dy = d/dv, one has exactly ^(r) = qJ 2 . Consequently, given 
that the coefficient q = completely determines all scalar linear iter¬ 
ative equations in normal form and with second-order source equation 
fl2.9l) . it follows from Corollary [1] that for any order n and under the 
substitution fl4.12l) . one must have 

'I'LrIyI = y‘”' + + ■ ■ • + /ijy = o, (4.i3a) 

where the scalar equation 

yi-) + Aly^--^'> + --- + Aly = b, (4.13b) 

is the corresponding nth order iterative equation with source equation 
02.91) . and more precisely one has = ^d^iq]. That is, the coefficients 
Ai^ in the iterative vector equations of the form 04.13a|) are obtained 
explicitly as differential functions of q. In the sequel diag(xi,... ,Xp) 
will denote a diagonal matrix with diagonal entries xi^... ,Xp. 

Although the matrices R in 04.lip have been calculated explicitly 
only for the dimension m = 2, it is however easy to verify directly that 
for any square matrix of order m > 1 of the form 

R = diag {ul ,..., , or R = u^ {vij ), 

where u and the Uj for j = 1,... m are arbitrary nonzero solutions of the 
second-order source equation 02. 9 p while the rij are scalars, the matrix 
R commutes with its derivative and >I(r) = qJm. Consequently, for all 
orders n > 1 any iterative equation generated by the corresponding 
operator 4/nor of the form 04.8p lies in the canonical class. 

Due to the fact that the normal form of equations in the canonical 
class consists of uncoupled iterative equations, their general solutions 
can also be obtained from those of the second-order source equation 
02. 9 p by a simple superposition principle. Indeed, the general solution 
w = (wi,..., Wm) of a system of linear equations of the form 04.lap is 
given simply in terms of two linearly independent solutions u and v of 
<m by 

n 

Wi = , for i = 1,..., m, (4-14) 

i=i 



















12 


where the Cij are arbitrary scalars. Indeed, this follows simply from 
the fact that fid.lap consists of identical scalar equations of the form 
fid.lbp , the general solution of which is given by fl4.14p for a fixed value 
of h In fact, as demonstrated in [22] it suffices to know a single nonzero 
solution u of the source equation fl2.9l) to derive the general solution of 
fl4.1ap . as the general solution w can indeed also be expressed as 


Wi = 


/ , CijU 

i=i 


n—1 


dx V 
^ ) 


for i = 1,..., m. 


(4.15) 


Similarly, for a general equation of the form (I3.2p in the canonical 
class, once we know the matrix Q in (14. 5 p dehning the transformation 
mapping it to its normal form, the general solution of (13. 2 p is simply 
given by y = Qw, with w given by (I4.14p . 


5. Concluding remarks 

We have shown in this paper that any system of linear equations re¬ 
ducible by an invertible point transformation to the canonical form is 
iterative, and determined arbitrary classes of iteration operators gener¬ 
ating equations in the canonical class. However, we’ve also noted that 
contrary to the scalar case proved in [16], not all iterative systems of 
linear equations belong to the canonical class. It is therefore an open 
problem to hnd the most general form of operators 'h = R^ -|- S whose 
generated equations all lie in the canonical class. 

If it was known as it is in the case of scalar equations that a system 
of linear equations is reducible by an invertible point transformation to 
the canonical equation if and only if it is of maximal symmetry, then the 
result of this paper asserting that any linear equation in the canonical 
class is iterative would be an extension of that obtained by Krause 
and Michel [16] for scalar equations. However, to our knowledge it has 
been established only for linear systems of orders two [1] and three 
[2] that the canonical class and the class of linear systems of maximal 
symmetry are the same. 
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